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We study the effect of topology for a random matrix model of QCD at nonzero imaginary chemical 
potential or nonzero temperature. Nonuniversal fluctuations of Dirac eigenvalues lead to normal- 
ization factors that contribute to the 6 dependence of the partition function. These normalization 
factors have to be canceled in order to reproduce the 9 dependence of the QCD partition function. 
The reason for this behavior is that the topological domain of the Dirac spectrum (the region of the 
Dirac spectrum that is sensitive to the topological charge) extends beyond the microscopic domain 
at nonzero imaginary chemical potential or temperature. Such behavior could persist in certain 
lattice formulations of QCD. 



I. INTRODUCTION 

Chiral random matrix models pj (2] have been used 
with great success in the past 15 years to understand and 
compute universal features of the QCD Dirac spectrum; 
see Ref. for a review. The effect of temperature or 
chemical potential can be included in these models in a 
schematic way to obtain qualitative, nonuniversal results 
for the QCD phase diagram. The main purpose of this 
paper is to point out and clarify a number of subtleties 
and nonuniversal effects that can arise when the effects 
of topology become important in such schematic random 
matrix models. In particular, we shall see that nontrivial 
normalization factors can occur which will be related to 
nonuniversal properties of the Dirac spectrum. 

Let us first consider QCD at zero temperature with Nf 
quark flavors, which for simplicity we take to be of equal 
mass m. The QCD partition function, Z^'~"^ , can be 
considered at fixed 9 angle or at fixed topological charge 
ly. In the former case, the 6 angle can be introduced 
according to (see, e.g., 
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where ma (tol) is the mass that couples right-handed 
(left-handed) quarks with antiquarks of opposite chiral- 
ity. We assume m to be real and positive. 

If the number of right-handed and left-handed modes 
differs by i^, the product of the fermion determinants re- 
sults in an overall factor e*''®, and we have 



(2) 



This relation can be inverted to give the QCD partition 
function at fixed v, 
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which corresponds to a path integral restricted to gauge 
fields of topological charge ly. 



It is generally assumed that the gauge field mea- 
sure does not depend on the topological charge. When 
topological excitations can be considered as independent 
events, the central limit theorem dictates that the distri- 
bution of topological charge is given by 
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where V is the space-time volume and x is the topo- 
logical susceptibility at = 0. In the quenched theory 
(or, equivalently, for heavy quarks), X = Xg is a mass- 
independent constant, whereas for light quarks, the topo- 
logical charge is screened, resulting in a topological sus- 
ceptibility at = given by [5] 



(5) 



where S is the absolute value of the chiral condensate for 
m = and = 0. 

In the microscopic domain of QCD, where the Comp- 
ton wavelength of the pion is much larger than the size of 
the box, the mass and Q dependence of the QCD partition 
function is given by a random matrix theory (RMT) with 
the same global symmetries as those of QCD. Contrary to 
QCD, random matrix partition functions are defined in 
terms of integrals over the matrix elements of the Dirac 
operator at fixed topological charge rather than integrals 
over gauge fields at fixed angle, which contain the sum 
over topological charges. In this paper we will study ran- 
dom matrix theories that are deformed by an imaginary 
chemical potential or temperature. The deformation pa- 
rameter will be denoted by u. 

Given a random matrix partition function at fixed z/, 
the partition function at fixed d is defined by 
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where we separated a potentially nontrivial normaliza- 
tion factor Mu and a weight factor Vv from Z^'^'^(to). 
The factor Vy corresponds to the quenched distribution 
of topological charge given in Eq. (|4| (with x = Xq)- 
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Other {/-dependent normalization factors that may arise 
in random matrix models of the QCD partition function 
are included in the factor My. One of our objectives is 
to discuss the significance of these two factors. We shall 
see in Sec. Ill A| that, contrary to QCD or chiral random 
matrix theories at m = [U [6] , AC can become a nontriv- 
ial function of the deformation parameter. On the other 
hand, as will be shown in Sec. |III B| for light quarks it 
makes no difference whether or not Vu is included in the 
sum over v [6J. 

A related question we would like to address in this 
paper is which part of the Dirac spectrum is sensitive 
to the topological charge. The answer to this question 
could depend on the parameters of QCD or the chiral 
random matrix model, and we shall see below that this 
is actually the case. It also depends on the value of 
the quark mass, for which we distinguish the following 
scales: (i) The microscopic scale [H, where mVH is 
kept fixed in the thermodynamic limit. This corresponds 
to the e regime of chiral perturbation theory [7 . (ii) The 
chiral scale where m^/V is kept fixed in the thermody- 
namic limit. This corresponds to the p regime [8^ of chiral 
perturbation theory, (iii) The macroscopic domain with 
m ~ Aqcd- In the microscopic domain, the mass depen- 
dence of the QCD partition function is given by chiral 
random matrix theory. Actually, this domain extends 
beyond the microscopic domain all the way to the chiral 
scale. Therefore, it is appropriate to borrow the name 
"ergodic domain" from the theory of disordered systems 
[9] to distinguish the domain m <C 1/AqcdV^ from the 
microscopic scaling domain. Note that we will sometimes 
consider the limit where mVT. approaches infinity with 
the understanding that the thermodynamic limit is taken 
first so that m is still in the microscopic domain. 

The issues that will be addressed in this paper are al- 
ready manifest for one quark fiavor, and for simplicity 
we will only discuss this case. The one-fiavor QCD parti- 
tion function, given by the average fermion determinant, 
is a function of the quark mass and of the 9 angle or 
the topological charge v. If the eigenvalues of the (anti- 
Hermitian) Dirac operator at fixed v are denoted by lA^, 
the QCD partition function at fixed v can be expressed 
as 



(7) 



where the average is over gauge fields with fixed v. 

We know that in the microscopic domain (and in fact 
in the ergodic domain) the mass dependence of the one- 
flavor QCD partition function in the sector of topological 
charge is given by HI |7] 



Z,(m) - /.(mFS) 
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For large values of the argument the modified Bessel 
function I^, becomes insensitive to its index v, and thus 
Eq. (|8| implies that the average fermion determinant 
does not depend on the topological charge when mVS 3> 



1 . In terms of Dirac eigenvalues one way to realize this is 
when only eigenvalues below this mass scale are affected 
by topology [see Eq. (|7|]. However, more exotic scenar- 
ios are also possible. It could be that eigenvalues beyond 
the microscopic domain are sensitive to the topological 
charge. If m is in the microscopic domain, this might 
result in a i^-dependent overall factor that could de- 
pend on the deformation parameter u and restores the v- 
independence of Zi, for mVH ^ 1. To find out whether 
this scenario is realized, it makes sense to introduce the 
notion of the topological domain of the Dirac spectrum, 
which we define to be the part of the Dirac spectrum that 
is sensitive to the topological charge. 

In QCD we have A^ = 1 and, from Eqs. ^ and (|8|, 
the universal 6 dependence of the partition function is 
given by 
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It is plausible that the standard scenario discussed after 
Eq. Q applies in this case, i.e., the topological domain 
of the Dirac spectrum does not extend beyond the micro- 
scopic domain. Exotic scenarios such as the one discussed 
above could occur in certain lattice formulations of QCD, 
and it would be interesting to test this directly. We shall 
further comment on this point in the conclusions. 

The ergodic domain of QCD is given by random ma- 
trix theory, but since the average fermion determinant is 
sensitive to all eigenvalues, it could be that deformations 
of the random matrix model result in a topological do- 
main that extends beyond the microscopic domain. In 
this paper we will see that this may happen in random 
matrix models at nonzero temperature/imaginary chem- 
ical potential. 

The 6 dependence of random matrix theories at non- 
zero temperature was discussed before in the literature 
[lOJ. In that work the temperature was introduced such 
that it only affects the eigenmodes corresponding to non- 
zero Dirac eigenvalues. This resulted in the same 9 
dependence as in the zero-temperature random matrix 
model [1 . Among others it was shown that the 0{m?) 
term in the chiral Ward identity does not contribute in 
the chiral limit. This is not always the case. It was 
recently shown in the framework of chiral perturbation 
theory that in the superfiuid phase of QCD at nonzero 
chemical potentials the 0{m?) term cannot be neglected 
prj . In this paper we will see that the 0{m?) term in 
the chiral Ward identity contributes to the topological 
susceptibility for random matrix partition functions at 
nonzero temperature/imaginary chemical potential if the 
M-dependent normalization factor AO is not included. 

The structure of this paper is as follows. Chiral ran- 
dom matrix theories at zero and nonzero deformation 
parameter will be introduced in Sec. [IT] The random ma- 
trix models are solved in Sec. |III| where we also discuss 
the normalization factor A^ and the distribution of the 
topological charge Vu- In Sec. |IV| we show that the chiral 
condensate for one flavor only has the correct behavior 
if the normalization factor A^ is included. The origin of 
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J\fy is studied in Sees. [V| and |VI[ In Sec.|V]we show that 
it is related to the extent of the topological domain, and 
in Sec. |VI| we find that the contribution of the pseudo- 
scalar susceptibility does not vanish if Mu is not included. 



Concluding remarks are made in Sec. VII 



II. CHIRAL RANDOM MATRIX MODELS 



A. Definition of tlie random matrix model 



The random matrix model for Nj = 1 in the sector of 



topological charge v is defined by [IJ 
Z^{m) = Cn,u j VWdei{D + m) e" 



with the random matrix Dirac operator defined by 

iW^ 



D = 



iW^ 



(10) 



(11) 



The integral "DW is over the real and imaginary parts of 
the elements of the random matrix W , which has dimen- 



sion p X g. The Dirac operator (111 has \p — q\ exact zero 
modes. For this reason we interpret 

u^p-q (12) 

as the topological charge. The total number of modes 

N = p+q (13) 

will be interpreted as the volume. This corresponds to 
the choice of mode density 
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The normalization factor Cj^^^, is chosen such that the 
quenched partition function is normalized to unity, i.e.. 
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We will consider this random matrix model in the pres- 
ence of an imaginary chemical potential iu. Using the 
chiral representation of the 7 matrices, the u-deformed 
Dirac operator is given by [TH |T3J [T3] |TS] 



D{u) 





iW^ + m 1 



iW + iu 1^ 




(16) 



where (Ipxg)^^ — ^ke- Alternatively, u can be interpreted 
as a schematic temperature as was done in |T5J [T31 [H] . 
The argument goes as follows. The temperature enters 
in the Dirac operator through the matrix elements cor- 
responding to Oq, with eigenvalues that are given by the 
Matsubara frequencies. We include only the temperature 



dependence given by the lowest two Matsubara frequen- 
cies by adding the p x q temperature matrix iT to iW 
and iW^ in Eq. (11 1, where 



for k < min{p, q}/2 , 
for k > min{p, g}/2, 



(17) 



and — for k £. Using the invariance of the 
integration measure under unitary transformations W 
UWV~'^ with U e U(p) and V € U(g), the temperature 
matrix can be transformed into a diagonal matrix with 
all diagonal matrix elements equal to u, so that the Dirac 
operator is given by Eq. (16 1. 



In the following, we shall refer to the model defined by 
Eq. ( 16 ) as model A. 



B. Other random matrix models 



Equation (16 1 is not the only way to introduce a non- 



zero temperature. Another possibility [TU] is to first par- 
tition the N modes into Nq = p + q "zero" modes and a 
fixed number A^i of "nonzero" modes, with |^| = |p — 9] 
actual zero modes of the Dirac operator. An A^i x Ni tem- 
perature matrix is then added to the nonzero-mode com- 
ponent of the Dirac operator, while the zero-mode matrix 
elements remain temperature independent. In terms of 
the Dirac operator ( [T6| this means that we add to an 
{Ni/2 + p) X {Ni/2 + q) random matrix W a diagonal 
matrix with A'^i /2 elements equal to iu and min{p, q} el- 
ements equal to zero. (This is technically equivalent to 
the model considered in Ref. [T^, although the physics 
background is different.) In the following, we shall refer 
to this model as model B. 

A third possibility is to add to W a, random matrix 
with matrix elements that are proportional to u. This 
model was introduced in Ref. (16j for imaginary u (i.e., 
real chemical potential) to describe the microscopic do- 
main of QCD at nonzero baryon chemical potential. For 
real u, this results in a model that differs from the orig- 
inal model (111 simply by a rescaling of the parameter 
E according to S ^ S/Vl + u^. This model will be re- 
ferred to as model C. Note that this model does not have 
a chiral phase transition. A less trivial model is obtained 
by introducing two or more different imaginary chemical 
potentials [17], but we will not discuss this possibility in 
this paper. 



III. SOLUTION OF THE RANDOM MATRIX 
MODELS AND NORMALIZATION FACTORS 



In this section we solve the random matrix models that 
were introduced in the previous section. We will find that 
the universal 9 dependence is not recovered for model 
A at u =/= unless additional normalization factors are 
included. 
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Solution of model A 



we find the universal result [H [B] 



In this subsection we solve the random matrix model 



A given by Eq. (10 1 with Dirac operator (16 1. The pro- 



cedure is standardjsee, e.g., [DII^]). We start by writing 
the determinant as a Grassmann integral, 



det(i:'(u) + m) 



d^pd-ip exp 



(18) 



and perform the Gaussian average over the random ma- 
trix elements. After a Hubbard-Stratonovich transfor- 
mation and integration over the Grassmann variables we 
obtain the following a model: 



dada* (^l + u^la + m\ ^ 
X (a + m)P{a* + m)« e-(i/2)JVS=a^* ^ (^g) 



where n — min{p, g}. Notice that the ^/-dependent nor- 
malization constant introduced in Eq. (10) has canceled. 
After changing variables a —^ a — m and a* —^ a* — m 



in Eq. ( 19 1 and then expressing the integral over {a, a*) 
in polar coordinates (r, if) , the angular integral results in 
a modified Bessel function, and the partition function is 
given by the remaining integral over r, 

Z^{m) = 27r Hdr I,{mN^\)r\''\+\r^ + y2)(JV-|H)/2 

^0 

X e-(i/2)A^s^(r^W) (20) 

For large A^, this partition function can be evaluated 
by a saddle-point approximation. For m in the er- 
godic domain, the saddle point in the broken phase is 
at = 1 /S^ — u^. To leading order in 1 /N the partition 
function is given by 

Z^'-^(m)^/,(mAS^(w))Tl=^l, (21) 

where irrelevant prefactors have been ignored and 

S^(w) = St(u) with t{u) = Vl - S2u2 . (22) 

A second-order phase transition to the chirally symmetric 
phase occurs a,t Uc — [T^ . 

The 9 dependence of the partition function is obtained 
after performing the sum over v according to Eq. (|6|. 
As will be explained in detail in the next subsection, for 
light quarks the sum is not affected by the distribution 
function Vi, [S]. We will therefore set ~ I. 

Let us first consider the case u = and take Afi, = 1. 
Using the identity for Bessel functions given by f |18j. 
Eq. (9.6.33)) 



E 



(23) 
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This shows that we do not need nontrivial normalization 
factors at u — Q. 

Now consider the case u 7^ 0. Because of the factor r''^! , 
in this case Eq. ( 21 1 depends on v for miVS^(M) > 1. 



This is a nonuniversal result and would also lead to a 
nonuniversal 9 dependence of after summing over v. 
However, these problems can be fixed by introducing a 
M-dependent normalization factor 



(25) 



Then with the replacement E — s- T,^ (u) the sum over v is 
the same as for u — 0. Again the sum is not affected by 
the distribution function Vu, and we find the universal 
result 



(26) 



In QCD an imaginary chemical potential is equivalent to 
a constant vector field and can be gauged into the tempo- 
ral boundary conditions of the fermion fields. This is not 
the case in random matrix theory, and therefore it should 
not come as a surprise that we need a :^-dependent nor- 
malization factor to recover the correct 9 dependence. In 
agreement with universality properties of Dirac spectra 
at fixed v [T7J [ini [Ml HI] this normalization factor does 
not depend on the quark mass. 

When u approaches Uc = 1/S, higher-order terms in 
the saddle-point approximation of Eq. ( 20 ) become im- 
portant, and the integral has to be performed exactly. 
We will not further elaborate on this and only discuss the 
parameter domain where the leading-order saddle-point 
approximation is appropriate. 

We will discuss further properties of model A in later 
sections but first turn to a discussion of the necessity of 
Vi, and to a comparison with models B and C, where no 
M-dependent normalization factors will be needed. 



B. On the necessity of Vi^ 

For large at fixed x the modified Bessel function 
can be approximated by ([18], Eq. (9.3.1)) 



ix/2) 



(27) 



Therefore, if m is in the microscopic domain, the sum 
over v in Eq. ^ is convergent without the Gaussian fac- 
tor ^. 

The sum over v can be performed, up to exponentially 
suppressed contributions, using the approximation ^ 



(28) 
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which follows from the uniform large-order expansion of 
the modified Bessel function and is valid for 1 ^ <C a; 
f|18j. Eq. (9.7.7)). It makes no difference whether or not 
we include the factor Vv in Eq. (|6| since 



,-(i/V2W)((l/™S(u)) + (l/x,)) 



^-{u^ /2mNT.{u)) 



(29) 



for m in the ergodic domain. The topological suscepti- 
bility at 6* = is therefore given by Eq. ([s]). From the 
approximation ( [28| we also see that all topological sectors 
with <^ mNYXu) contribute equally to the partition 
function. 

It was argued by Damgaard 6 that the factor Vi, 
should be absent in the sum over v in Eq. (|6]), although 
he also pointed out that the quenched limit could not be 
taken properly in this case. Our point of view is that 
the presence of Vu is immaterial for m in the microscopic 
domain, but that Vu becomes important at length scales 
below the inverse rj' mass where it is believed to deter- 
mine the local topological susceptibility and leads to the 
Witten-Veneziano formula for the rf mass [251E51l2^l25j . 
Beyond this scale the topological susceptibility at = 
is given by Eq. 



C. Comparison with models B and C 

For fixed topological charge v the partition function of 
model B is given by 

Z^{m) = j dGda*{\a + mf + u^)^^'"^ 

X (cr + m)P(cr* + m)9e-(i/2)Ars=<T<T* ^ 
or, after introducing polar coordinates, 



Z^(m) = 27r / dr Iu(rnNY?r)r^'>"^^ (r^ + v?) 
Jo 

,-(l/2)AfS^(r^+m^) 



2\Ni/2 



X e 



(31) 



Note that this partition function becomes independent 
of 1^ for large mNTi. Since the correct 9 dependence is 
obtained at m = this model does not require additional 
normalization factors. The sum over v with V^, = 1 re- 
sults in 



Z^(m,^)=27r / dre™^^ '■™^^^°+i(r2 + t.2^ 
Jo 

-(l/2)JVS^(r^+m^) 



2\Nl/2 



X e 



(32) 



Using a saddle-point approximation for large N, we find 
the universal 6 dependence 



^mNi:"{u) cose 



where il3| 



1 - S2u2 + ^(l + E2^2)2_4J]2y2^^/^' 



(33) 
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(34) 



The partition function of model C at deformation pa- 
rameter u is equivalent to Eq. ( 19 1 at u = with 
E 5]'-^(m) = E/Vl-l-w^, and we thus have the uni- 
versal result 



(35) 



Hence we see that in order to obtain the universal 9 
dependence of the partition function neither model B nor 
model C requires normalization factors AO that depend 
on the deformation parameter u. However, let us em- 
phasize again that the appearance of the AO is a generic 
feature in RMT. In the remainder of this paper we will 
identify mechanisms that are responsible for this behav- 



IV. CHIRAL CONDENSATE AND TOPOLOGY 

The case Nf = 1 we address in this paper is special 
since there is no SU(A^/) x SU(7V/) symmetry that could 
be spontaneously broken. Nevertheless, there could still 
be a chiral condensate, which can be calculated in the 
usual way, 



|(V'V')| = :^9mlogZ(m,( 



(36) 



The parameter S introduced earlier is defined to be equal 
to I (V'V') I at 6* = for TO ^ and V oo. The functions 
E(u) computed in Sees 
w-dependent chiral cone 



III A and III C correspond to the 

These 



: ensate in the same limits 
limits can be taken in different orders pS], either 

lim lim -a„JogZ(TO, 61 = 0) (37) 
or in the reverse order 

S(2) = hm hm ^d,n log Z{m, 9 = 0). (38) 



In Eq. ( 37 1 , a nonzero chiral condensate implies the 



breaking of the Uyi(l) symmetry by instantons or the 
chiral anomaly [271 , whereas in Eq. ( 38 1 a nonzero chiral 



condensate implies "spontaneous symmetry breaking" in 
the following sense. At fixed topology the QCD partition 
function has a Ua(1) symmetry (in fact a covariance ex- 
cept at v = where we have a symmetry). A nonzero 
chiral condensate spontaneously breaks this U/i(l) sym- 
metry at fixed topology. 

From the universal expression ^ for the one-fiavor 
partition function it is clear that the order of limits 
should not matter. We will now see that for model A 
this is only the case if the normalization factors AO are 
included. Because in this section we only consider model 
A we omit the superscript A. Using Eq. ^ and the mass 
dependence of Zu{m) given by Eq. (21 1, we find that 'E^^'> 
of model A is given by 



E*-"'"-' = lim lim 



d,n[J\fiZi{m)+J\f^iZ^i{m)] 



, (39) 
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where the factor Vu has dropped out of numerator and 
denominator since it is essentially constant for v <C ^fN . 



Using the result (21 1, we obtain 

= (M/AAo)I]r2 = (AAi/AAo)E(l - ll'^u') . (40) 

Next we calculate the chiral condensate using the re- 
verse order of limits. Bas ed on t he discussion in Sec. lIIIBl 
we find that for |zy| ^ \/mNY, the condensate for fixed 
V does not depend on v. Its value is therefore equal to 
the value in the v — Q sector. This was calculated in 
Ref. T2|, resulting in 
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(41) 
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We thus see that the two condensates are only equal 
if the normalization factor Ni /Aq = 1/ \/l ~ S^u^ from 
Eq. ( |25| is included. Therefore we have a nice consistency 
check of Eq. (p5l. 



So far, we have assumed that we can choose p and q to 
have arbitrary v = p — q. Let us now fix the total number 
of states N . In this case the Dirac operator with v zero 
modes has nonzero off-diagonal blocks of dimension (A^+ 
v)/2 X {N - i/)/2; see Eqs. ^ and (fTsl. This imphes 



that the parity of the topology is the same as the parity of 
N. In the following we assume that N, and therefore also 
I/, is even. Equation ^ shows that the chiral condensate 
can be extracted from 



(E(i))^ 



lim lim ^^^^-^ 
N^oom.^oN^ Z{m,e = 0) 



0) 



(42) 



For m — > 0, the numerator receives contributions from 
1/ = and 1/ = ±2, while only the v ^ sector con- 
tributes to the denominator. For the ly — contribution 
we find |4] 



(S(i))2^o = lim lim^ 



1 dlZoim) 



IN/2 
\k=l 



1™ AT9 \ 



(43) 



FIG. 1: Convergence of the ratio R„ of determinants for i/ = 
2 and = as a function of the number n of eigenvalues 
included for an ensemble of 10® 400 x 400 matrices. Results 
are shown for u = 0.0, u — 0.5, and u — 0.8. The shaded 
areas correspond to the statistical errors. 



( 43 1 and ( 44 1 evaluate to 



1. 



.(1)n2 _ -^^2_4 



2 

l„o AAf2 



(45) 
(46) 



so that the correct result for the chiral condensate is ob- 
tained with the normalization factors from Eq. ( p5| . 

The question we wish to address in the sections below 
is why model A requires the u-dependent normalization 
factors ( 25 1 . We will relate this question to the properties 
of the Dirac eigenvalues. As we have already discussed 
in the introduction, the requirement that Z^(m) be inde- 
pendent of v for mVTi ^ 1 can explain why a normaliza- 
tion factor Ml, ^ 1 is needed if the topological domain of 
the Dirac spectrum extends beyond the microscopic do- 
main. Equation (44) shows that the consistency relation 
E^^) = E^^) should also be related to the properties of 
the Dirac eigenvalues, to which we turn now. 



Nf = l 



where the average includes the fermion determinant. The 
right-hand side of Eq. ( 43 1 is dominated by the smallest 
eigenvalues. Note that this contribution is independent 
of the normalization of the partition function. The con- 
tributions of = ±2 to the condensate are the same and 
can be written in terms of the Dirac eigenvalues as 



=a\2 



(44) 



where averages without subscript are with respect to the 
quenched partition function. This is essentially the ratio 
of the fermion determinants in the sectors v = 2 and 
1/ = 0. In the random matrix model A the expressions 



V. EIGENVALUE FLUCTUATIONS AND 
MICROSCOPIC UNIVERSALITY 

In the numerical calculation of this section we keep 
N fixed as discussed at the end of the previous section. 
Motivated by Eq. ( 44 ) , we consider the ratio i?„ of the 
products of eigenvalues for v = 2 and v = Q as a. function 
of the number of eigenvalues included in the product. 



Rn = 



1 



iV2 



nLi(Ar°)^ 



(47) 



For n = iV/2 all eigenvalues are included in the prod- 
uct, and for model A the value of this ratio follows from 
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FIG. 2; Topological shift AA„ of the eigenvalues for an en- 
semble of lO'' 400 X 400 matrices (top) and an ensemble of 10^ 
800 X 800 matrices (bottom). The shaded areas correspond 
to the statistical errors. 



Eqs. (441 and (46 1 



N^oa ' 8 



(48) 



We have evaluated the ratio i?„ numerically for model 
A, using an ensemble of 10^ random matrices (16) of 



dimension N — 400 distributed according to the Gaus- 
sian factor in Eq. ( [To| . The mass has been set to zero. 
In Fig. [l]we plot the ratio Rn/Roo versus n for u = 0, 
u — 0.5, and u — 0.8. We observe that for u = the 
ratio of determinants saturates in the ergodic domain 
^ VN = 20). This is not the case for u = 0.5 and 
u = 0.8, where all eigenvalues contribute to the ratio of 
the two determinants. 

This is further illustrated in Fig. [2] where we plot the 
ratio 



{K~") ^ (K+i) 



(49) 



versus n. The motivation for constructing this particu- 
lar ratio is as follows. The microscopic eigenvalues are 
expected to behave universally after rescaling with the 
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FIG. 3: Convergence of the ratio Rn of determinants for u = 2 
and = as a function of the number n of eigenvalues in- 
cluded. Results are for an ensemble of 10^ 400 x 400 matrices. 
The shaded areas correspond to the statistical errors. 



chiral condensate and the volume. The universal result 
for the spectral density of microscopic eigenvalues in the 
quenched case and in the topological sector i' is [28] 



PsiO = I [MO' - J.+iiOJu-iiO] , (50) 



where J^y is a Bessel function and ^ = XNY.. Its large-f 
behavior is given by 



1 cos(i/7r - 2^) 
PslO = ^ 



(51) 



so that for A„A^S > 1 we have (K^^) ~ (K+i) and 
therefore AA„ 0. This is indeed what we find in Fig. [2] 
for u — 0. Notice that Eq. ( [50| has been obtained by 
taking the microscopic limit and is only valid for eigen- 
values well below the chiral scale. For w = we find 
that AA„ = also beyond the microscopic domain and 
conclude that in this case the topological domain does 
not extend beyond the microscopic domain. For u =/= 0, 
however, the situation is completely different. All eigen- 
values are in the topological domain and only the first 
few eigenvalues show universal behavior. Comparing the 
results for N = 400 and N = 800 in Fig. (S) we observe 
that the universal domain, i.e., the domain where the 
eigenvalue ratio AA„ does not depend on u, increases 
with N proportional to \/N. This is in agreement with 
microscopic universality for u < Uc = 1/S, which states 
that the distribution of low-lying eigenvalues is universal 
after rescaling them by the chiral condensate. If we con- 
sider the Dirac spectrum around x, the correction terms 
to this universal behavior are of the order Nx^. This im- 
plies that the number of eigenvalues with universal fluc- 
tuations around A = scales with \/N. 

Based on Fig. [2] a plausible explanation for the behav- 
ior of the ratio of the determinants seen in Fig. [l]can be 
given in terms of the u dependence of the average posi- 
tion of the eigenvalues. For this reason we plot in Fig. [3] 
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FIG. 4: Number variance versus the average number n of 
levels in the interval [0, n] for an ensemble of 10® 400 x 400 
matrices. The curves for u — 0.0 and u = 0.5 only start to 
deviate from each other at n > 15. 



FIG. 5: Topological shift AA„ and the ratio Rn/Roo for an 
ensemble of 10*^ 400 x 400 matrices with Ni/N = 0.75 for 
model B [TD]. 



the same ratios as in Fig. [T] but normalized with respect 
to the average positions of the eigenvalues. The ratio i?„ 
defined by 

- (nr! (ArV(Ar))') 

Rn = A (52) 

(nLi(Ar7(Ar°)) ) 

is shown for u = 0.0, u = 0.5, and u — 0.8. 

We conclude that the u dependence of the ratio of the 
determinants is almost exclusively due to the effect of u 
on the average position of the eigenvalues. 

In the theory of disordered systems, a frequently used 
measure to test the breakdown of universality is the num- 
ber variance [9J. This is the variance of the number of 
levels in an interval containing h eigenvalues on average. 
In Fig. [4] we display the number variance versus the 
average number n of eigenvalues in an interval starting 
at zero. The curves for u = 0.0 and u = 0.5 coincide for 
n < 15, while the deviations between u = 0.0 and u — 0.8 
are already significant for n > 5. This is in agreement 
with the discussion of Fig. [2j 

In Fig. [5] we show the behavior of the Dirac eigenvalues 



in model B. We observe that in this model the topological 
domain does not extend beyond the microscopic domain 
even for m ^ 0. This is also the case for model C, which 
at u ^ is equivalent to model A at u = after rescaling 
the chiral condensate S T,^{u). The results for model 
C are therefore identical to the u — results in Figs.[T]-|4] 
We thus have a further piece of evidence that nontrivial 
normalization factors AO only appear if the topological 
domain extends beyond the microscopic domain. 



VI. TOPOLOGICAL AND PSEUDOSCALAR 
SUSCEPTIBILITY 

As mentioned in the introduction, the 9 dependence of 
the QCD partition function is obtained by introducing 
left-handed and right-handed quark masses according to 
z — me*^ and z* = me~^^ , respectively; see Eq. ([T]). De- 
noting the left-hand side of Eq. Q by Z{z, z*), with the 
superscript QCD omitted for simplicity, the topological 
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susceptibility at arbitrary angle is given by 
Xi = \^{{^^)-{^?)=-\^dl\ogZ{z,z*) 



{zd, + z*d,,)logZ{z,z*) 



(53) 



1 

V 



z'^dl + z*'^dl, ~ 2zz*d,d,> 



log Z{z,z*) 



Because mdm — zd^ + z*dz', the first term on the 
right-hand side of this equation is equal to m| <^'0^/'^ | ; 
see Eq. (36). The second term on the right-hand side 



of Eq. (53 1 is equal to m? times the pseudoscalar (PS) 
susceptibility given by 



- V{zi!L^R - Z*4!RljjLfj^_^ . 



(54) 



(55) 



Thus Eq. ( 53 1 becomes 

Xt = m|(-0V')| + nn?Xvs ■ 

This is the well-known chiral Ward identity relating xt to 
the chiral condensate and the pseudoscalar susceptibihty 
[5]. Note that = 5^ cos 61 0(m). 

The random matrix partition function Z"^(to, 0) with 
TO in the ergodic domain can be calculated explicitly from 
Eq. MjSetting = 1 according to the discussion in 
Sec. '" " ■ 



IIIB 



We will set Mv = t~^'^^'^^~^\ where setting e 
to or 1 allows us to switch between including or not 
including N^- 



We first replace the Bessel function in Eq. ( 20 1 by 
the integral representation 



1 



ii/ip-\-x cos ip 



(56) 



sum the resulting geometric series in v, and perform a 
saddle-point approximation of the radial integral includ- 
ing next-to-leading order corrections in m to find 



Z^(to, 



27r 



dtp 



1 



r2£ 



/o 1 — cos ip + t 
X exp [mN'ET cos{ip — d)] 
1 



2e 



X exp 



4r2 



Note that 



1 



r2£ 



u^o 27r(l — 2r'^ cos ip + t^^) 



but also 



r2£ 



lim — ^— 

£-•0 27r(l — 2t^ cos ip + t^^) 



(57) 



(58) 



(59) 



Therefore for e ^ or m 

Z^(TO,6')~exp niNT.TCOSt 



we find 
1 



4t2 



A^m^E^cos^6l 



(60) 



and thus by Eq. ( 53 I 

xt{u) =mY.^{u) cos e + 0{m^), (61) 

which is consistent with results obtained by Crewther [5] . 
We conclude that for u = or if we include the normal- 



ization factor ( 25 ) the contribution of the pseudoscalar 



susceptibility vanishes in the chiral limit. 

The situation is different, however, if we do not in- 
clude the Afp. For mNH 3> 1 the contribution of the 
pseudoscalar susceptibility to the topological susceptibil- 
ity becomes comparable to that of the chiral condensate 
but with opposite sign and thus leads to a significant sup- 
pression of the topological susceptibihty (see Fig. [6]). Be- 
cause the saddle-point approximation breaks down close 
to w = 1 we do not plot the curves of Fig. [6l in this re- 
gion. For mNH <C 1 the exponent in Eq. ( |57[ ) can be 
expanded, and after evaluating the integral analytically 
we find 



Z^(m, 9)^1 + mNJ:T^+' cos 9 . 



(62) 



This result agrees with Fig. |6] and shows that in this 
limit the contribution of the pseudoscalar susceptibility 
at u 7^ is small also without Afi, . 

Metlitski and Zhitnitsky have recently found another 
situation in which the 0{m?) term in Eq. (55 1 becomes 
important, i.e., the superfluid phase of QCD with two or 
three colors JJJ. Of course, if we include the AC in model 
A (as we should) we do not see this effect. Nevertheless, 
our observation may potentially be of importance; see 
the conclusions. 

For models B and C no normalization factors AO are 
needed to ensure a vanishing contribution of the pseudo- 
scalar susceptibility. 

The vanishing of the contribution of the pseudo- 
scalar susceptibility also imposes constraints on the v- 
dependence of pseudoscalar correlators and can be used 
as a check of results that were recently derived for the e 
domain [2^150]. 



VII. CONCLUSIONS 

It is well-known that random matrix models for QCD 
at zero imaginary chemical potential (or temperature) u 
have the correct 9 dependence. In this paper we have 
shown that this is not automatically the case for m ^ 0. 
We obtain the correct 9 dependence only after introduc- 
ing i/-dependent normalization factors AC in the sum over 
topologies. 

To explain this we have introduced the topological do- 
main of the Dirac spectrum, which is defined as the part 
of the Dirac spectrum that is sensitive to the topological 
charge. We have shown that for tt = the topological do- 
main coincides with the microscopic domain. This is also 
the case at u for models for which no i/-dependent 
normalization factors are needed to obtain the correct 9 
dependence. However, for the model we analyzed that 
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FIG. 6: Contributions to for ^ = 10^ and m = 5 • 10 ^ (top) and m = 5 ■ 10 ^ (bottom) with or without the normahzation 
factor A/L. We set S = and E = 1. The curves were obtained by numerical evaluation of Eq. ([ST]) in connection with Eq. ([53|. 



requires nontrivial normalization factors, the complete 
Dirac spectrum is inside the topological domain. This 
results in a partition function that gives universal be- 
havior for small Dirac eigenvalues, but has bulk spectral 
correlations that depend both on u and on the topologi- 
cal charge. In the thermodynamic limit this leads to an 
additional u-dependent factor in the partition function 
at fixed topological charge which results in an incorrect 
6 dependence of the partition function. To obtain a par- 
tition function with the usual behavior in the chiral limit, 
one has to introduce additional j^-dependent normaliza- 
tion factors in the sum over topologies. 

Our observations are of potential importance for lattice 
QCD at nonzero imaginary chemical potential or tem- 
perature. Depending on, e.g., the fermion formulation or 
the algorithm used, it could be that nontrivial normal- 



ization factors are needed in the sum over topological 
sectors, and these could even persist in the continuum 
limit. To find out whether such normalization factors 
might be necessary, it would be interesting to determine 
the topological domain as a function of the deformation 
parameters. This is feasible with current lattice tech- 
nology. To be consistent with the general properties of 
QCD, the topological domain should not extend beyond 
the microscopic domain. Future work will tell us if this 
interesting picture prevails. 
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